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ABSTRACT 



We utilize high-resolution (R~60,000), high S/N (~100) spectroscopy of 17 
cool Pleiades dwarfs to examine the confounding star-to-star scatter in the A6707 
Li I line strengths in this young cluster. Our Pleiads, selected for their small pro- 
jected rotational velocity and modest chromospheric emission, evince substantial 
scatter in the linestrengths of A6707 Li I feature that is absent in the A7699 K I 
resonance line. The Li I scatter is not correlated with that in the high-excitation 
A7774 O I feature, and the magnitude of the former is greater than the latter 
despite the larger temperature sensitivity of the O I feature. These results sug- 
gest that systematic errors in linestrength measurements due to blending, color 
(or color-based T eS ) errors, or line formation effects related to an overlying chro- 
mosphere are not the principal source of Li I scatter in our stars. There do 
exist analytic spot models that can produce, via line formation effects, the ob- 
served Li scatter without introducing scatter in the K I line strengths or the 
color-magnitude diagram. However, these models predict factor of >3 differ- 
ences in abundances derived from the subordinate A6104 and resonance A6707 
Li I features; we find no difference in the abundances determined from these two 
features. These analytic spot models also predict CN line strengths significantly 
larger than we observe in our spectra. The simplest explanation of the Li, K, 
CN, and photometric data is that there must be a real abundance component to 
the Pleiades Li dispersion. We suggest that this real abundance component is 
the manifestation of relic differences in erstwhile pre-main-sequence Li burning 
caused by effects of surface activity on stellar structure. We discuss observational 
predictions of these effects, which may be related to other anomalous stellar phe- 
nomena. 
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Introduction 



Dramatic differences in the Li abundances of main sequence stars in open clusters 
stand in stark contrast to the greater uniformity that is the general rule for many other 



elem ents. T 
early iHerbig 



lecomp 



(119651 ) 



exity of the observed pattern of st ellar Li depletion was recognized 



Wallerstein. Herbig fc Contil ( 119651 ) and can be traced to the fragility 



of the species. Lithium is destroyed by proton capture at relatively low stellar interior 
temperatures (of order 2.6 million K for typical densities); these conditions are achieved 
for most low mass stars during the pre-main sequence (pre-MS) phase, which yields a 
predicted mass-dependent depletion pattern (e.g., Iben 1965). Physical processes negle cted 



1965|). 



in standard stellar models can also induce lithium depletion ( IWeymann fc Searsl 
Open cluster studies have revealed two other generic features of Li abundance patterns: 
the existence of a dispersion in abundance at fixed mass, composition, and age and the 
existence of main sequence depletion even in stars with convection zones too shallow to be 
able to burn lithium (for reviews see Pinsonneault 1997 and Jeffries 2006). 



The specific case of star-to-star Li dispersion in the Pleiades, and other young c 



such as IC 2602 ( IRandich et al. 



20011 ) and a Per (IBalachandran. Lambert Stauffer 



usters 



1996|), 



has been extremely challenging to understand from both theoretical and observational 
perspectives. A substantial dispersio n in the equivalen t widt h of the 6707 A Li I resonance 



line for cool Pleia ds was reported 



larger data set of 



Soderblom et al. 



subordinate Li I feature by 



( 



Duncan fc Jonesl (119831 ) and confirmed by the much 



1993). Subsequ e nt ob servations of the weak 6104 A 



Ford. Jeffries k, Smalleyl (120021 ) yielded a consistent result. The 



star-to-star Li equivalent width variations, which are superposed on a strong T e g-dependent 
depletion pattern that is presumably governed by pre-MS Li destruction in the deep surface 
convection zones whose extent is determined uniquely by stellar mass for a given composition 
and age in standard stellar models, have been difficult to explain in the context of differing 
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surface Li abundance s. Star-to-star Li abundance variations can develop in the context 



of rotational mixing (IPinsonneault. Kawaler fc Demarque 



and are certainly seen in older open clu sters such as M67 (IPasquini. Randich &: Pallavicini 



990 



Ryan fc Del ivannis 



1995) 



1997 



Jones. Fischer &: Soderblom 



19991 ) . However, little such mixing is expected to have 
occurred at the Pleiades age (~100 Myr). Furthermore, the highest rate of mixing is 
expected in rapid rotators, while the rapid rotators tend to populate the upper envelope of 
the Li equivalent width distribution in the Pleiades. 

It thus remains important to establish whether the star-to-star range in Pleiades' Li 
equivalent widths reflect a real star-to-star range in Li surface abundances before strong 
conclusions can b e drawn about the und erlying physical cause of the Li line strength 



dispersion. While 



Soderblom et al. 



(119931 ) suggest the Pleiades Li dispersion reflects real 



abundance differences, they note the difficulty for this explan ation that com es from the 
scatter they find in the Pleiades' A7699 K I resonance feature. iJeffriea (Il999l ) confirm this 
K I scatter, and note the concomittant danger in concludi ng that the Pleiades Li d ispers ion 



ste ms from real ab u ndan ce variations-a caveat echoed by IStuik. Bruls fc Rutted (119971 ) 



and 



line formation, 
corre 



Carlsson et al.l (|1994l). who note t he similarity in the details of the Li and K resonance 
King fc Schulerl (120041 ) catalog a rich variety of Li-K and alkali-activity 



King. Krishnamurthi fc Pinsonneault 



ations in young clusters, seen in the Pleiades by 
(120001 ). associating some 90% of the variance in Li and K in M34 (200-250 Myr) and IC 
2391 (50 Myr) with the spread in chromospheric emission (though this emission is likely 
a proxy for surface inhomogeneities in the ca se of M34). Because rapid rotators in the 



Pleia des evince apparent Li overabundances (ISoderblom et al. 



1993 



Garcia Lopez et al. 



19941 ). chromospheric emission may simply serve as a 



proxy for rotation in t his case. In 



Margheim et al 



(120021 ) find that 



their study of the 100 Myr old NGC 2451 A and B, 
significant Li scatter, reminiscent of that in the Pleiades, is spuriously produced as an 
artifact of deriving Li abundances from equivalent widths; they suggest that accounting (via 
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spectral synthesis) for blending features in rotationally broadened spectra might eliminate 
the Pleiades Li scatter. 

Here, we seek to address the important question of a real Li dispersion component in 
cool Pleiads by utilizing new spectroscopy of higher resolution and/or S/N than in most 
previous studies to examine Li and K spreads in slow projected rotators that demonstrate 
modest chromosperic emission relative to other cluster members. We conclude that there is 
evidence for a real dispersion in abundance at fixed effective temperature, and (in the final 
section) advance the idea that variations in pre-main sequence depletion stemming from 
differences in stellar physical parameters arising from surface inhomogeneities, rather than 
variations in the rate of mixing, may be implicated. 



2. Observational Data and Analysis 

We selected 17 Pleiads having a range of T e g believed to evince significant Li dispersion 
(Soderblom et al. 1993; King, Krishnamurthi & Pinsonneault 2000, hereafter KKP) and 
having v sin i<12 km s -1 . All objects are radial velocity members and have at least one 
proper-motion study indicating cluster memb ership. Fourteen of the objects show no radial 



19921). H II 571 and 2406 



velocity or photometric evidence of binarity (jMermilliod et al. 

are single-lined binaries, but were included because their mass function determination 

and photometric decomposition suggest companion y-band contamination at only the 

<1% level. H II 298 is the lower mass 6-7 arcsec distant visual binary companion to H II 

299. Our sample is listed in Table 1, which gives information on S/N, luminosity, color, 

membership, and rotational velocity. Figure 1 shows our objects in the Ca II near-IR 

triplet- and Ha-based chromospheric emission versus V a planes defined by the Pleiades 

sample from KKP. Our projected slow rotators exhibit well-below maximal chromospheric 

emission values at a given luminosity. Tab. 1 

Fig. 1 
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We used the Hobby-Eberly Telescope 9.2-m telescope and its High Resolution 
Spectrograph (HRS) to obtain spectroscopy of our targets on numerous nights from August 
2002 to October 2003. Wavelength coverage from 5095 to 8860 A was achieved over the 
HRS 2-CCD mosaic. The S/N in the Li I A6707 region is 80-160 per pixel. The 0.5 arcsec 
slit width yielded a nominal resolving power of i?~60,000. Standard reductions were 
carried out with the I RAF package to accomplish bias removal, scattered light subtraction, 
flat-fielding, order extraction, and wavelength calibration. We do not conduct spectrum 
synthesis of the A6707 Li I and A7699 K I lines since quantitative deblending is not required 
for these features in our spectra and absolute abundances are not the principal topic of 
interest. Instead, we focus on equivalent width differences of these lines for any subset of 
the program stars having a narrow range of color. 

Equivalent widths of the A6707 Li I and A7699 K I resonance lines were mea- 
sured with Guassian and V oigt profile fitting routines using the SPECTRE package 



( IFitzpatrick &: Sneden 



19871 ). The principal potential danger in this approach is the 
mild blending some 0.4A to the blue of the Li I feature. However, its presence is of 
little importance here: the small projected rotational velocities mean that any residual 
contamination in the empirically deblended Li equivalent widths given in Table 1 is limited 
to a few ml (the size of random uncertainties in the line measurement), and the strength 
of the blending Fe I (augmented by CN in cooler dwarfs) feature proves to be nearly 
invariant among Pleiads of given intrinsic color. More importantly, even the full blend 
contribution is fractionally small, achieving strengths of 15-20 mA only in our coolest stars 
whose Li equivalent widths are an order of magnitude larger. A conservative upper limit on 
equivalent width uncertainties including continuum placement is <10 mA. 

For the purpose of full disclosure, we note that two spectra of H II 152 taken 14 
October 2003 UT are clearly at least double-lined. Two other spectra (08 December 2002 
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UT) are single-lined. The flux levels and gross spectral appearance do not clearly establish 
a (queue) target misidentification for the 2003 spectra. However, this seems the simplest 
explanation inasmuch as several previous independent high-resolution spectroscopic studies 
(e.g., Soderblom et al. 1993; Wilden et al. 2002; Boesgaard, Armengaud fc King 2003) 



have n ot noted an SB2 classification. While the BY Dra classification by 



Kholopov et al 



( 119891 ) is not inconsistent with binarity, it appears to be based on the modest photometr ic 



variation attributed to 4. Id variations from spot pattern migration by 



Magniitskvi (119871 ). 



We proceed here with the 2002 spectra and the assumption that the spectrum is single-lined. 



3. Results and Discussion 



Our results are presented in Figure 2, containing purely observational planes showing 
equivalent widths versus (deredenned) (B — V) color. Several conclusions can be reached: 
First, our new data nearly remove the star-to-star scatter in the A7699 K I equivalent 



widths seen in 



Soderblom et al. 



therefore important in the 



(1993 ) for our sta rs. While measurement uncertainty is 



Soderblom et al. 



(119931 ) results, this can not be the source of 



most of the star-to-star Li dispersion in their data since this scatter (up to a factor of 2 in 
equivalent width for 0.72<(£> — V) <0.84) persists in our own higher- quality data. Second, 



the significant Li scat ter in our projected s 



ow rotators indicates that abundance errors due 



20021 ) do not provide the primary explanation of the 



to blending features ( IMargheim et al 
Pleiades Li scatter. 

Finally, the vast extent of scatter in Li I compared to K I suggests that simple color 
or T e ff errors or the influence of an overlying chromosphere are not a significant source of 
star-to-star Li scatter in our Pleiads. Color or temperature errors would lead to a similarly 
large dispersion in the similarly temperature-dependent K I lines that is absent in our 
data. Since the details of A6707 Li I and A7699 line formation in cool dwarfs are similar 



Fig. 2 
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( IHoudebine fc Doyle 



1995 



Stuik. Bruls fc Rutten 



19971 ). the dramatic difference between 



the Li and K in our Pleiads rules out differences in the global properties of overlying 
chromospheres as a sole or dominant source of the Li scatter. Indeed, our sample selection 
was made to mitigate such effects: as seen in Figure 1, our stars evince modest scatter 
in chromospheric emission compared to a more representative cluster sample. This is not 
to say that such differ ences are unimportant contributors to Li scatter in other samples 



( IKing &: Schuler 



20041 ) . Indeed, even here the K I equivalent width of our reddest Pleiad 



(H II 298), which also possesses the largest chromospheric emission index in our sample, is 
anomalously large compared to a modest linear extrapolation o f the K I- ( B-V) relation that 
can be seen to be linear as red as (B — 1/)q~1.4 in Figure 1 of iRandichl (120011 ) for the 35 
Myr IC 2602 cluster. 



3.1. Photospheric Inhomogeneities 



It is expected that surface magnetic activity (spots and plages) could alter Li I line 
strengths. Spatially resolved solar observations show variations by fa ctors of ~2 and t ens 



in the Li I equivalen t width in spot and plage regions, respectively ( iGiampapa 



Patterer et al. 



19841 ). 



(119931 ) have found Li I line strength variations in weak-lined T Tau stars 



that are not correlated with chromospheric emission vari ations, but a re con s istent in size 



with t hose expected from the simple spot/plage model of iGiampapal (119841 ) 



Jeffries et al 



(119941 ) find rotational modulation of the A6707 (and/or nearby blended features) in the 



young Local Association K6 dwarf BD+22 4409. 

One might think that the stark difference between the factor of ~2 scatter in A6707 
Li I line strengths for 0.7<(B — y)o<0.8 versus the absence of scatter in A7699 K I over 
this same range is a powerful argument against spots as a source of Li scatter in our 
stars. Surprisingly, this is not necessarily the case. The influence of various analytic spot 
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models on the location of 
seen in Figures 3 and 6 of 



leiads in the K,Li line str e ngth versus color plane can be 



Barrado y Navascues et al 



(1200 ll ). The models with the two 



lowest photosphere-spot temperature contrasts (their models 2a and 2b) and 80% spot 
coverage are able to displace stars redward and to higher Li I equivalent width into the 
0.7<(B — V)o<0.8 color range in such a way as to introduce a near factor of 2 scatter in Li 
line strength, but move stars nearly parallel to the intrinsic K I line strength versus color 
locus. I.e., the analytic spot models can, in fact, produce the Li I scatter we observe while 
not introducing substantial scatter in the K line strengths. 



As seen in Figure 2 of iBarrado y Navascues et al.l (120011 ). these models also move the 
Pleiads roughly parallel to the main-sequence in the H-R diagram such that significant 
photometric scatter, AV>0.1 mag, is not introduced either; such scatter is not evinced 
by our stars (Figure 3). Barrado y Navascues et al.'s (2001) larger photosphere-spot 
temperature contrast models (models 2c and 2d) can reproduce the observed scatter in Li 
line strengths with spot coverage fractions of ~60%; however, their Figure 2 shows that 
such stars would be subluminous by AV^>1 mag in the color- magnitude diagram. The 
photometric data in Figure 3 clearly excludes these large temperature contrast spot models. 
The large photosphere-spot temperature contrast models only modestly alter the star's 
color (Figure 2 of Barrado et al. 2001). As a result, a similar degree of scatter is introduced 
into the Li and K linestrengths (Figures 3 and 8 of Barrado et al. 2001)-a prediction also 
excluded by our spectroscopic data. 

In sum, the marked difference between K I and Li I scatter can not exclude a spot 
origin for the Li scatter in our Pleiads. While the indistinguishable spread of Li in 
the (V — I)- and V^-Li EW planes (Figure 4) compared to that in the (B — F)-based 
plane of Figure 2 does not betray the presence of spots , two spectroscopic signatu r es do 



robustly exclude the analytic spot models 2a and 2b of 



Barrado y Navascues et al 



(1200 lh 



Fig. 3 
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as a source of the Li scatter. The first is the difference in Li abundances derived from 
the A6707 Li I resonance line and the weaker blend ed 1.85 eV A6103.6 Li I subord inate 
feature(s) previously investigated in the Pleiades by Ford. Jeffries fc Smalley ( 2OO2L We 



selected several interesting sets of Pleaids of similar color but with significantly different 
A6707 Li I line strengths (Table 2). Effective temperatures an d microturbulent velocitie s 



are ad apted from the spectroscopic parameters determined by iFord. Jeffries fc Smalley 



(120021 ). Following these authors, we assumed log g = 4.5 for all stars since the derived Li 



abundances are insensitive to the assumed gravity. We determined Li abundances from the 



A6707 resonance line 
program discussed in 



rom our measured equivalent widths via Dr. A. Steinhauer's LI FIND 



King fc Schulerl (120051 ); these are given in column 6 of Table 2. 



Abundances from the A6104 Li I features were determined vi a spectral syn thesis 



carried out with an updated version of the LTE analysis MOOG 



Snedenl (119731 ) and 



Kurucz model atmospheresj. The A6104 region linelist was formed utilizing atomic data 

I, 



from the Vien na Atomic Line Databa se (IKupka et al. 



CN data from iDavis fc Phillips 



2000) and Kurucz line datazl, and 



( 119631 ) . The linelist was calibrated by adjusting some 



oscillator stren gths in order to produce solar syntheses matching t 



re Kurucz solar flux atlas 



( Kurucz 



20051 ) using input abundances from lAnders fc Grevessd (Il989l ) except for CNO 



solar values of log N(C)= 8. 3 9, log N( N)= 7.78, and log iV(0)= 8.68 were a d opted from 



Allende Prieto. Lambert fc Asplund 



(I200lh . The 



Asplund. Grevesse fc Sauvall (120051 ) and 
Pleiades syntheses assumed input s caled solar abundances with [m/H]= —0.06 based on the 



results of 



Ford. Jeffries fc Smallevl (J200J) 



Sample A6104 Li I spectra and syntheses are shown in Figure 5. Abundances are given 
in the penultimate column of Table 2, whose final column lists the A6104- and A6707-based 



Fig. 4 
Tab. 2 



x http : / /kurucz . cf a. harvard.edu / grids, html 
2 http:/ /kurucz. cfa.harvard.edu/linelists. html 
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Li abundance differences. There are two notable features in the abundance comparisons. 
First, the significant Li scatter persists when measured from the A6104 feature. Second, 
we find no significant difference between the resonance line- and subordinate line-based 
abundances. The mean difference (6104—6707) is +0.01 dex with small scatter (±0.08 dex, 
s.d.). Fig. 5 



We find that the 



Barrado y Navascues et al. 



(1200 ll ) analytic spot models reproducing 



the observed Li dispersion of our sample in the range 0.7<(B — ^)o<0.8, lead to significantly 



larger predicted 6104—6707 L i differences 



latter using the linelist from 



King et al 



Synt 



reses of the A6104 and A6707 Li I regions-the 



( 119971 ) updated with VALD atomic data and 



recalibrated to the solar flux spectrum- were performed using the photosphere and spot 



param eters of models 2a and 2b with 80% spot coverage from 



Barrado y Navascues et al 



( 120011 ) assuming log g = 4.5, £ = 2.0 km s" 1 , [m/H]= —0.06, and an input abundance of log 
iV(Li)= 2.80. The photospheric and spot spectra were weighted by their respective Planck 
functions and coverage fractions before adding and renormalizing. 

The resulting A6104 and A6707 photosphere+spot spectra were analyzed as observed 
data via spectrum synthesis computed for T e g- = 5225 (the value of the spotted models 
consistent with flux conservation), log g = 4.5, £ = 2.0, and [m/H]= —0.06. The A6707 Li 
abundances deduced from the spotted synthetic spectrum are log iV(Li)= 2.90 and 3.02 for 
spot models 2a and 2b, respectively. The corresponding A6104 Li abundances are 3.32 and 
3.62. The Li differences (6104—6707 A) deduced from the two Li I lines- 0.42 and 0.6 dex 
for spot models 2a and 2b-are significantly larger than the zero difference exhibited by our 
stars. 

The second spectroscopic signature that is inconsistent with the analytic spot models 
reproducing the observed Li spread in our stars is the strength of CN lines in both the 
A6104 and A6707 region. The CN line depths in the T e s = 5225 spotted synthetic spectra 
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are some 3-5 times deeper than observed in our actual object spectra. Figure 6 illustrates 
a related important point: while the warmer photosphere (T p h „ t = 5655 K; 20% coverage 



Barrado y Navascues et al 



and co ol spots (T spot = 4870 K; 80% coverage) in model 2a of 
(120011 ) conspire to yield a spotted star of T e ff~5225, the spectrum of such a model is not 
equivalent to that of an unspotted model with T e g = 5225. Figure 6 indicates the CN 
features are significantly stronger in the spotted model. 



Fig. 6 



3.2. Over-Excitation and -Ionization Effects 



Schuler et al. 



(120041 ) find that our cool Pleiads' O abundances derived from the 



high-excitation O I A7774 triplet in our spectra show a d ramatic increase with declining T e g 



and significant star-to-star scatter. 



Schuler et al. 



(j2006al ) show that photospheric hot spots 



provide a plausible explanation for this behavior. In order to examine whether the Li and O 
scatter are related, the A6707-based Li and A7774-based O abundances were fit as functions 



of T e fj using low order polynomials. To ensure consistency for this purpose 



Schuler et al 



i abu ndances 
(12004); these 



Schuler et al 



(I2004h O 



were determined with LIFIND using the stellar parameters of 
were updated for H II 298 using our (B — V)o value, adjusting the 
abundance (to [0/H]= +0.80) in the process. The resulting abundance residuals (observed 
minus fitted), shown in Figure 7, are not correlated. Moreover, the rms dispersion of the 
observed Li abundances (~0.22 dex) is twice that of the O abundances (~0.12 dex) despite 
the fact that the temperature sensitivity of the (logarithmic) O abundances is 2-3 times 
larger than that for Li. We conclude that the dispersion in our Pleiads' Li is not associated 
with that in O I. Fie. 7 
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4. Conclusions, a Proposed Explanation, and Future Work 

4.1. Summary and Key Conclusions 

We have utilized high-resolution and -S/N spectroscopy to measure A6707 Li I and 
A7699 K I line strengths and A6707- and A6104-based Li abundances in 17 slowly (projected) 
rotating cool Pleiades dwarfs. A significant factor of ~2 dispersion is seen in the A6707 
Li I line strengths over the 0.72<(I? — V^)o<0.82 color range; this large Li scatter is also 
inferred from A6104 subordinate line-based abundances. The scatter in our selected sample 
eliminates line blending due to rapid projected rotation as a source of Li dispersion in 
our Pleiads. In contrast to previous studies, our high-resolution and S/N data evince no 
substantial scatter in the A7699 K I line strengths of our particular Pleiad sample. This 
stark distinction relative to the Li linestrengths excludes simple color or T e g errors or line 
formation effects due to an overlying chromosphere as the source of Li scatter in our stars. 

The difference in the dispersions of the K I and Li I line strengths does not, however, 
exclu de spots as a source of the Li scat ter; in particular, the analytic spot models 2a and 



2b of 



Barrado y Navascues et al. 



(120011 ) can produce the Li scatter we observe without 
leading to significant scatter in K I line strengths or in the color magnitude diagram. The 
equivalence of the A6707- and A6104-based Li abundances in our stars does, however, 
exclude these spot models, which predict a factor of >3 difference in the resonance 
and subordinate line-based Li abundances. These spot models also would lead to CN 
linestrengths significantly larger than observed in our spectra. 

The simplest explanation of the Li I dispersion in our Pleiads is scatter due to real 
abundance differences. There are numerous candidate mechanism(s) to explain these 
pre- or near-ZAMS abundance differences: chemical inhomogenieties, magnetic fielc 



differences, variable mass accretion, and individual rotational histories (jVentura et al. 



1998 
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Garcia Lopez et al.lll994f ). While the effects of chromospheres on Li I line 
effects of r otation on measuremen ts of the Li I line strength espoused by 



or mat ion and the 



f l2004h and 



Margheim et al 



King fc Schuler 



( 120021 ) may contribute to an (illusory) Li dispersion, our results 
suggest that there must also be a real component to the Pleiades Li dispersion-at least over 
the range 0.7<(B — V)q<0.8 seen in Figure 2. 

The size of the real component is difficult to determine directly; however, future 
work can deduce the size and mass range of real scatter since sources of possible illusory 
dispersion are amenable to observation. Systematic effects of rotation on Li abundance 
measurement of rapid rotators or from low spectral resolution (neither characterize our 
data) can be mitigated by determining abundances from spectrum syntheses with suitable 
accounting of macroscopic broadening. The influence of chromospheric emission and surface 
magnetic activity (spots and plages) on Li I (and K I) lines can be measured by searching 
for correlated temporal variations in alkali line strength and chromospheric emission and 
photometric ind ices via a sim ultaneous spectroscopic and photometric monitoring progra m 



of the Pleiades ((Jeffries 



19991 ). Our Li equivalent widths and those of ISoderblom et al. 



(I1993[ ) for 5 stars (H II 152, 263, 2126, 2311, and 2366) differ at the >2a level, providing 
impetus for such a program. Additional constraints may be gleaned from observations of 
9 Be in our cool Pleiads, as well as in similarly cool Hy ades stars, whose degre e of Li scatter 



is unknown due to vanishingly weak Li line strengths (ISoderblom et al. 



19951 ) 



4.2. A Proposed Differential Pre-MS Li Burning Mechanism 

We believe that the most likely theoretical explanation of a real star-to-star Li 
dispersion in the Pleiades is a range in pre-MS lithium depletion. Such a range would 
naturally emerge if the radii of protostars of the same mass, composition, and age during the 
epoch of pre-MS lithium depletion were somehow different. To explain the empirical trend, 



- 16 - 



one would also require the more rapidly rotating and active stars to have been larger. There 
is now emerging evidence th at activity impac t s the radii of both main sequence and pre-MS 



stars in precisely this sense. 



Torres &: Ribad ( 120021 ) found that the radii of the near-twin 



stars in YY Gem were far too large to be consistent with the predictions of standa rd stellar 



Morales. Ribas 



Jordil (J2008) found a 



Berger et al. 



(120061 ) 



theory. Subsequent data confirmed this result, and 
correlation between activity and the radius excess in low mass stars, 
also found this to be a relatively common phenomenon in their measurements of the radii 
of field M stars, although it is more challe nging; to perform a rigorous test in the absence 



of direct mass constraints in such systems. 



Andronov fc Pinsonneault 



(120041 ) . in a study of 



cataclysmic variable systems, found that excess activity could induce radius changes of the 
proper order of magnitude (see also Chabrier, Gallardo, & Baraffe 2007). 

The physical mechanism linking activity-related radii differences with those in Li 
depletion is that spots inhibit convective energy transport, requiring the star to carry 
the flux through a smaller effective volume. The star therefore requires a larger radius, 
corresponding to an effective reduction in the mixing length or efficiency of energy transport. 
This in turn reduces pre-MS lithium burning because the central temperature and pressure 
of a star is decreased if the radius is increased. We performed a simple numerical test of this 
phenomenon. A 10% difference in radius (of order seen in active eclipsing binaries) during 
standard pre-MS convective burning for a 0.9 M solar abundance model led to a factor 
of 2.5 increase in the remaining surface lithium abundance at an age of 100 Myr relative 
to an inactive model with the same mass but a smaller radius. How the spot properties 
of single pre-MS stars (the precursors of our Pleiad sample) compare to those of active 
binaries remains an open question. Some meaningful context, however, might be provided 
by chromospheric emission fluxes. We note that the log i?'(HK) indices for active solar-type 
binaries widely range from —4.4 to —3.3 at a given color or Rossby number (Figure 4 of 
Montes et al. 1996), while the corresponding Pleiades (pseudo-) indices (which must be 
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transformed from Ha and Ca II infrared triplet data) widely range from -4.5 to -3.5 (Figure 
3 of King et al. 2003); i.e., active binaries show wide chromospheric emission at the same 
levels and with similar wide dispersion as present-day Pleiads. 

When/if spot filling factors are reduced subsequent to pre-MS Li depletion, stars 
would converge to similar radii, but would retain a fossil record of their differences in 
earlier stages. If more spotted stars on the main sequence retained their lower effective 
temperatures, one would also shift more massive (and less depleted) stars to the same 
effective temperature as less massive (and more depleted) bare stars, further amplifying the 
appar ent Li dispersion. There is ev i dence that activity can differentially impact protostellar 



radii. IStassun. Mathieu fc Valentil ( 120071 ) found a reversal in the predicted T e ff-radius 
relationship between a more active primary and a less active secondary of a brown dwarf 
binary system in the young (1 Myr) Orion Nebula Cluster. 

Our suggested explanation of the Pleiades Li dispersion has several attractive features. 
It explains why the dispersion is large amongst cooler stars, which experience significant 
pre-MS depletion, compared to hotter stars, and the propensity for more rapidly rotating 
(and perhaps more spotted) cool dwarfs to exhibit larger Li line strengths. It also can 
qualitatively explain the apparent narrowing of dispersion in cool main sequence stars 
in older open clusters as follows. Beyond the Pleiades age, more heavily spotted and 
rapdily rotating stars deplete more Li than barer and more slowly rotating stars due to 
rot at ionally- induced main-sequence mixing; as spottedness declines with increasing age 
and stellar radii concomitantly relax, the former stars also become hotter relative to the 
latter stars. The combined effect would be a reduction in the star-to-star dispersion in cool 
cluster Li abundances that is controlled by the timescales of stellar angular momentum 
loss and spot coverage dimunition. Because the former timescale could, in principle, be 
longer than the latter, the passage of additional time would then witness the resurgence 
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of substantial Li dispersion due to rot at ionally- induced main-sequence mixing acting in 
stars whose radii and temperatures are no longer scattered due to surface inhomogeneities. 
This is not inconsistent with the observation that the ratio of rapid and slow rotators 
in the ~100 Myr Pleiades is the inverse of that of Li- rich a nd -poor solar mass stars in 



the 4.6 Gyr cluster M 67 (jJones. Fischer fc Soderblom 



19991 ). Our suggested explanation, 



which recognizes the importance of surface inhomogeneities on stellar radii and possible 
differences in time evolution of rot at ionally- induced mixing and spot coverage, dovetail s 



Jones. Fischer fc Soderblom! (1199 



with the important observational picture painted by 
a significant Li dispersion in cool dwarfs in very young clusters such as the Pleiades that 
markedly declines in intermediate age clusters such as M34 and the Hyades, and then 
reappears in older clusters such as M 67. We caution, though, that the dwarf Li abundances 
in M 67 are available only for stars of T e g>5500 — 5600 K; abundances in lower mass 
M67 dwarfs akin to our cooler Pleiads have not yet been determined. It should also be 
noted that a variety of observational studies (e.g., Piau, Randich & Palla 2003; Randich 
et al. 2007, Jeffries et al. 2002) have suggested other mechanisms besides or in addition to 
rot at ionally- induced mixing as the source of main-sequence Li depletion and explaining the 



resulting main-sequence Li-T e ff morphologies manifested by open cluster data. 



4.3. Related Stellar Phenomena and Future Work 



It is also tempting to speculate that the starspot-radius mechanism might be connected 
with other puzzling stellar phenomena. The alleged e xistence o f a ga p (or gaps) in 



the number distribution of P o p I m ain sequence stars 



Mendozal (I1956h : 



Bohm-Vitense 



(119701 1; iRachford fc Canternal (120001 ) around (B — V)~0.3 could be accomodated if this 
color corresponds to the onset of surface spot formation driving stars to lower T e g values 
compared to bluer barer stars. The accompanying increase in stellar radii might give rise 
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to the surprising number of luminosity class 



region of the Hipparcos-based H-R diagram (INewberg fe Yanny 



I-IV stars that reside in the main-sequence 



19981 ). The precipitous 



rise in apparent [Fe II/H] abundances and more gradual but never theless surprising rise 
in apparent [Fe I/H] abundance in the Hyades with declining T e g ( lYong et al.l 120041 ). the 
rise in high-excitation permited [O I/H] abun dances with declining T e g in the Hyades and 



the even steeper rise in the younger Pleiades (ISchuler et al. 



line-based [O I/H] abundances with declining T e g in the 



Hyades (|Schi. 



2006a[). the rise in forbid den 



er et al. 



2006a), the 



20031 ) are at least 



rise in [Si I/H] abundances with declining T e g in M 34 (ISchuler et al. 
qualitatively consistent with a T e g- dependent disparity between spot-adjusted radii and 
assumed standard stellar radii. Such a disparity would lower log g with declining T e g; 
the sensitivity of the features noted above is such that lowering log g increases their line 
strengths in cool dwarfs. 

Future observational work is needed to confirm these speculative connections; 
quantifying or parameterizing spot coverage and properties rather than chromospheric 
emission indices per se that are used as measures of "activity" will be of particular 
importance. Such work includes: continued exploration of the reality of main-sequence gaps 
in open clusters, their evolution with age, and the spot properties of stars adjacent to these 
gaps; understanding the nature of luminosity class II-IV stars, in particular spot properties 
and physical radii, near the main-sequence region of the H-R diagram; the spot properties 
of stars as a function of T e $ in the Pleiades, M 34, and Hyades open clusters, and their 
association with the abundance anomalies (especially star-to-star variations in permitted 
line [O I/H] values in cool Pleiads and Hyads) listed above. 

A possible signature of the mechanism we propose here is the Li/Be ratio in stars 
sufficiently cool and old to suffer main-sequence mixing of 9 Be, but not so cool as to 
have suffered pre-MS convective burning of 9 Be. Such rot at ionally- induced main-sequence 
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mixing would deplete more Be (and Li) in stars with larger initial rotational velocities; 
however, if such stars were more heavily spotted, they suffered reduced pre-MS convective 
burning of Li; depending on the balance of spot-induced retarded pre-MS Li depletion 
and rotation-induced enhanced MS Li depletion, this effect could create a population of 
higher-Li/lower-Be stars whose number depends upon the initial distribution of rotation, 
its time evolution, and age. Identifying such stars using extant disk field abundance data is 
difficult given differences in stellar structure due to metallicity differences, unknown relative 
initial Li and Be abundances, and a dispersion in age. Such a search is best carried out in 
cool dwarfs of open clusters of at least intermediate (e.g., Hyades) age. For the purpose 
of constraining the origin of Li scatter in the Pleiades, a monitoring program to identify 
short- (days) and long- (years) term variations in Li I linestrengths and spot properties 
is critical. If our proposed origin of the Pleiades Li dispersion is correct, we expect to 
see variations in Li line strength that are consistent only with changes in photospheric 
structure/parameters and the details of spot-included radiative transfer alone; once these 
variations are accounted for, a significant dispersion in Li line strength should remain as a 
relic of spot-induced differences in pre-MS Li burning. While observationally challenging, 
these future observational programs (accompanied by detailed modeling of pre-MS Li 
depletion in stars of various spotted conditions) will be required to fully understand the 
magnitude, mass distribution, and thus the ultimate origin of real star-to-star Li dispersion 
in the Pleiades. 

JRK and SCS gratefully acknowledge support for this work from NSF awards 
AST-0086576 and AST-0239518, a generous grant from the Charles Curry Foundation to 
Clemson University, and a graduate scholarship award from the South Carolina Space 
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Fig. 1. — Our Pleiades sample (filled squares; Table 1) is plotted with other Pleiads from the 
study of KKP in the chromospheric emission (left: Ca II infrared triplet; right: Ha) versus 
V planes. The chromospheric fluxes are taken from Soderblom et al. (1993). 
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Fig. 2. — The top panels sh ow A6707 Li I (left) and A7699 K I (right) equivalent widths 



from Soderblom et al. 



(119931 ) versus dereddened (B — V) color. The bottom panels show the 
(empirically deblended in the case of Li) equivalent widths measured from our HET/HRS 
spectra. Representative error bars are shown in the upper left of each panel. 
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Fig. 3. — A color- magnitude diagram of our Pleiades stars. The > 1 mag scatter in V implied 



by lar ge photosphere-spot temperature contrast spot models of 



Barrado y Navascues et al. 



(120011 ) that might reproduce the Li line strength dispersion in Figure 2 is not present; any 
modest photometric scatter is not inconsistent with Barrado et al.'s smaller temperature 
contrast spot models, however. 
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Fig. 4. — The A6707 Li I equivalent widths of our Pleiades sample versus (V — I) color (left) 
and V magnitude (right). Typical uncertainties are shown in the upper left of both panels. 
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Fig. 5. — Observed (points) and synthetic (lines) spectra of the A6103.6 Li I region in H II 
250, 263, and 916 (top, middle, and bottom panels). The syntheses are shown for 3 input 
Li abundances for each star: no Li, the best fit Li abundance, and a Li abundance 0.3 dex 
larger than the best fit value. 
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Fig. 6. — (Dashed line) A sy nthetic spectrum of the A6707.8 Li I region for the 



Barrado y Navascues et al. 



(1200 ll ) analytic spot model 2a with T phot = 5655 K (20% cov- 
erage) and T spot = 4870 K (80% coverage), yielding a spotted model having T eff ~5225 K. 
(Solid line) A synthetic spectrum for T e g = 5225 K. All syntheses employ a 0.11 A FWHM 
gaussian smoothing and assume identical input abundances (solar) except for Li, for which 
log iV(Li)= 2.8 was utilized. The CN features are significantly stronger in the spotted model. 
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Fig. 7. — Residuals of Pleiades Li I abundances derived from our line strengths about a low- 
order Teg- dependent po lynomial versus the r esiduals in A7774-based O I abundances derived 
from the same spectra (ISchuler et al.ll2006al ). 
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Table 1. HET Pleiades Sample 



Star S/N V (B - V) (V - I) Mem. Probs. a RV b v sin i c EW(K I) EW(Li I) 



H II 


6707 










BFJ,SRSG,DH 


km s 1 


km s 1 


mA 




mA 


0152 


96 


10 


.60 


0.645 


0.685 


-0.98- 


5.2±0.2 


11.5 


194. 





191.2 


0193 


100 


11 


.18 


0.75 


0.80 


0.97,0.95,- 


7.6±0.2 


6.6 


216. 


4 


161.5 


0250 


97 


10 


.59 


0.645 


0.68 


-,0.98,0.57 


4.7±0.3 


6.4 


190. 


.6 


166.7 


0263 


133 


11 


.51 


0.84 


0.90 


0.98,0.98.0.17 


3.1±0.9 


7.8 


244. 


.7 


252.1 


0298 d 


80 


11 


.79 


0.89 


* 




4.8±0.2 


6.5 


355. 


.9 


266.8 


0571 


126 


11 


.14 


0.74 


0.85 


0.97,0.97,0.3 


5.7±0.1 


7.2 


223. 


,8 


164.0 


0746 


140 


11 


.18 


0.73 


0.82 


0.99,0.98,- 


6.5±0.4 


4.9 


208. 


.4 


107.6 


0916 


109 


11 


.59 


0.83 


0.89 


0.99,0.98,0.90 


4.2±0.7 


6.2 


241. 


3 


191.3 


1593 


128 


11 


.03 


0.73 


0.77 


0.99,0.04,0.01 


6.9±0.1 


2.4 


193. 


.5 


183.5 


2126 


119 


11 


.59 


0.81 
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0.98,0.00,0.05 
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244. 


.0 


184.9 


2284 
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.23 
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11 
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3 
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11 


.37 
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-0.98- 


6.2±0.2 


<5 


248. 


.4 


226.7 


2406 


150 


10 
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0.80 


0.99,0.96,- 


6.4±0.2 


5.2 


220. 


,8 


110.3 
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3179 
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.93 


0.53 
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a Proper motion-based membership probabilities from Ijonesl ( 119731 ). ISchilbach et al.l (119951 ). and 
Deacon fc Hamblvl (120041 ) 



b Mean radial velocities from WEBDA database: http:/ /obswww.unige.ch/webda 



^Rotational velocities (or conservative upper limits) from iQueloz et al.l ( 119981 ) 



d V and (B — V) from the photometric decomposition of 



Mermilliod et al. 



(1992) 
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Table 2. A6707 Li I versus A6104 Li I 



Star (B-V) T c $ £ EW(6707) A(Li) A(Li) AA(Li) 
HII K kins" 1 mA A6707 A6104 A6104 - A6707 



250 


0.645 


5726 


2.00 


166.7 


3.03 


3.09 


0.06 


74.fi 
t ^yj 


n 7i 


oooo 


9 1 n 


1 D7 (S 


9 ^1 




n D9 


1593 


0.73 


5383 


2.10 


183.5 


2.75 


2.73 


-0.02 


2284 


0.74 


5450 


2.88 


190.3 


2.87 


2.90 


0.03 


2311 


0.78 


5350 


2.30 


170.4 


2.64 


2.77 


0.13 


2366 


0.78 


5350 


2.30 


226.7 


2.98 


2.91 


-0.07 


2462 


0.80 


5283 


2.15 


110.3 


2.21 


2.37 


0.16 


2880 


0.82 


5250 


2.15 


142.3 


2.36 


2.30 


-0.06 


2126 


0.81 


5267 


2.15 


184.9 


2.61 


2.55 


-0.06 


916 


0.83 


5228 


2.15 


191.3 


2.60 


2.52 


-0.08 


263 


0.84 


5213 


2.13 


252.1 


2.97 


2.96 


-0.01 



